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NASA GLOBAL ATMOSPHERIC SAMPLING PROGRAM (GASP) 
DATA REPORT FOR TAPE VLOOO 1 

by James D. Holdeman and Erwin A. Lezberg 
Levis Research Center 


SUMMARY 

Atmospheric trace constituents in the upper troposphere 
and lover stratosphere are now being measured as part of the 
NASA Global Atmospheric Sampling Program (GASP), using fully 
automated air sampling systems on board commercial 747 
aircraft in routine airline service. Measurements of 
atmospheric ozone and related meteorological and flight 
information obtained during several GASP flights in March 
1975 are now available from the National Climatic Center, 
Asheville, North Carolina, In addition to the data from the 
aircraft, tropopause pressure data obtained from the 
National Meteorological Center (NHC) archives for the dates 
of the flights are included. This report is the first of a 
series of reports which will describe the data currently 
available from GASP, including flight routes and dates, 
instrumentation, the data processing procedure used, and 
data tape specifications. 


INTHODUCT ION 

This report announces the availability of ozone mixing 
ratio, static air temperature, and wind speed and direction 
data obtained at altitudes from 6 to 12 km during several 
flights of a Pan American World Airways 747 airliner 
(N655PA^ in March 1975. 

The objectives of the NASA Global Atmospheric Sampling 
Program are to provide baseline data of selected atmospheric 
constituents in the upper troposphere and lower stratosphere 
for the next 5-to-10 year period, and to document and 
analyze these data to assess potential adverse effects 
between aircraft exhaust emissions and the natural 
atmosphere. At present there is much uncertainty in 
environmental impact studies on this subject due to the lack 
of comprehensive, long-term upper atmospheric data (refs, 1 
and 2). 

The GASP program began in 1972 with a feasibility study 
ot the concept cf using commercial airliners in routine 
service to obtain atmospheric data. This program has 
progressed from design and acquisition of hardware (ref. 3) 
to collecting global data on a daily basis. Fully automated 



GASP systems are now operating on a United Airlines 747, a 
Pan American World Airways 747, and a Qantas Airways of 
Australia 747. The United airliner is collecting data over 
fhe contiguous United States and to Hawaii, Global coverage 
is provided by the Pan American and Qantas 747* s. Pan Am 
routes from the United States include a round-the-world in 
the Northern Hemisphere, trans- At lantic to Europe, 
trans-Pacif ic to the Far East, inter-continental to Central 
and South America, and occasionally trans-P acif ic to 
Australia. Hore frequent coverage in the Southern 
Hemisphere is provided by the Qantas 747 on transcontinental 
Australian flights and on flights from Australia to the 
South Pacific and Australia to Europe. The GASP system 
design, the measurement instruments, the on-board computer 
for automatic control and data management, and system 
maintenance procedures are described in reference 4. 

Although in the future GASP constituent monitoring will 
include ozone, water vapor, carbon monoxide, oxides of 
nitrogen, particle size and number density, halocarbons, and 
filter samples, the GASP data currently available are 
limited to ozone mixing ratios, and related meteorological 
and flight information from the aircraft system. Also 
included with the GASP data are tropooause pressure fields 
obtained from the National Meteorological Center (NHC) 
archives for the dates of the GASP flights. These data 
comprise GASP tape VL0001, which is available from the 
National Climatic Center, Asheville, NC, 23801. This report 
is the first in a series of reports which will present a 
catalog of data currently available from this program, Tn 
addition to announcing the data availability, this report 
describes the GASP instrumentation, data processing 
procedure, and data tape specifications. 


ROUTE STRUCTURE AND DATA ACQUISITION 

Flight routes foi. which data are given on GASP Tape 
VL0001 are shown on figure 1, All flights occurred between 
March 11 and March 31, 1975. The numbers tn parentheses on 
the figure indicate the number of flights for each route. 

On the tape, GASP data are grouped and identified by flights 
with the airports of departure and arrival designated by the 
standard three-letter airport codes (ref. 5). A listing of 
flights included in tape VL0CO1 by airport-pair, date, and 
data acquisition time, is qiven in table I. 

For each flight, data acquisition begins on asceht 
through the 6 km altitude flight level, and terminates on 
descent through 6 Km. A complete GASP sampling cycle is 60 
minutes, divided into 12 five minute segments. A 16 second 
recording is made at the end of each sampling segment. 

During alternate segments (at 10 minute intervals), air 
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sample data 
intervening 
calibration 


are recorded 
segments the 
modes to all 


instrument operation (if 
calibration mode is not 
instrument acquires air 


for ail instruments. During the 
system is in one of six different 
ow tor in-flight checks on 
required), Whenever any 
needed for a given instrument, that 
sample data during the segment. 


Cassette tapes, recorded in serial format, are removed 
from the aircraft at approximately tvo week intervals and 
transcribed to computer- compatible form for data reduction* 
At this stage, laboratory instrument calibration information 
required for data processing is included, redundant and 
non-usable data are removed, and the data are re-transcribed 
to final form and units. The detailed specifications and 
formats for the GASP data ate given in appendix A. Data for 
each flight begins with an FLHT record (table A-I) to 
provide flight identification information. This record is 
followed by a series of DATA records (table A-II), one for 
each recording made during the flight. 


MEASURE HE NTS 
Ozone 

Ozone measurements are made using a DASIBI Model 1003-AH 
continuous ultraviolet ozone photometer. The concentration 
of atmospheric ozone is determined by measuring the 
difference in intensity of an ultraviolet light beam which 
alternately passes through the sample gas and an ozone-free 
zero gas (generated within the instrument). The range of 
this instrument is from J to 20»000 ppbv (parts per billion 
by volume), with a sensitivity of 3 ppbv. This instrument 
is described in reference 6, and data from flight tests of 
the instrument are given in reference 7, The ozone 
instrument is checked at NASA-Lewis (over the range 0 to 
1000 ppbv) against an ozone generator which is calibrated by 
the one percent neutral buffeted potassium iodide (KI) 
method. The estimated accuracy of the KI procedure is seven 
percent. 

In-flight monitoring of the ozone i.istrument includes 
measurement of the instrument zero by flowing the sample 
through a charcoal filter external to the instrument, and 
measurement of the electronic span setting and control 
frequencies available from the instrument. The instrument 
is not calibrated in-flight with an ozone calibration gas 
due to the difficulty of generating a precisely known ozone 
concentration in the flight system. Periodic checks for 
calibration consistency are performed in the laboratory. 


The destruction of ozone in the Teflon sample lines from 
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the inlet probe to the instruoent, and in the Teflon-coated 
diaphragm pump that raises the sample pressure to 10 N/cm2 
(1 atn) , has been measured under conditions simulating 
operation in flight. Data from these tests can be expressed 
as 

03 = 03i ♦ 2.14 '/o3m (1 ) 

to within a standard deviation of 5.1 ppbv, where 03 is the 
ozone mixing ratio at the probe inlet, and 03m is the 
measured ozone mixing ratio. This relation has been used in 
the data reduction to correct the measured values to inlet 
(ambient) conditions. These data (03; see table A-II) are 
expressed in units of parts per billion by volume (ppbv). 


Plight Data 

In addition to the air sample measurements, aircraft 
flight data are obtained with each data recording to 
precisely describe conditions when the data are acquired. 
Aircraft position* heading, and the computed wind speed and 
direction are obtained from the inertial navigation system. 
Altitude, air speed, and static air temperature are 
collected from the central air data computer in the 
aircraft. Vertical acceleration information (an indication 
of turbulence) is taken from the aircraft flight recording 
system. Date and time are provided by a separate GASP 
clock-calendar unit. The formats and units for these data 
are given in table A-II, 


TROPOPAUSE PRESSURE DATA I 

The National Heteorological Center (NNC) is presently 
maintaining a library of gridded meteorological data fields 
accessible on various disk and magnetic tape systems (ref. 

8). Briefly, the data are interpolated to points on the N:.” 

65 X 65 grid, a square matrix map transformed from a polar 
steceographic map of the Northern Hemisphere, Among these 
gridded data are tropopause pressures, available on a twice 
daily basis (0000 and 1200 GHT) . 

Tropopause pressures are derived as a by product of the 
MHC objective analysis scheme which determines the height of 
constant pressure surfaces above each grid point. Vertical, 
mean layer temperature profiles, related directly to the 
vertical separation of the constant pressure levels, are 
calculated for each of the 4225 grid points, and fitted with 
a high order polynominal curve. By means of a slope testing 
routine, the tropopause is defined as the base of the lowest 
stable layer (pressures < 500 mb) within which the average 
lapse rate is < 2.5 degrees C/km. 
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The NMC tropopause pressure data arrays are included as a 
separate file (see appendix A) following the GASP data to 
provide indep.^ndent data tor analysis of GASP constituent 
neasureoents in terms of their tropospheric- stratospheric 
behavior. These TRPR records (table A-III) give the 
tropopause pressure at each of the 4225 coordinate points in 
the NBC matrix. The relations for obtaining the 
corresponding latitude and longitude from the NBC 
coordinates are given in appendix B. 


SELECTED ANALYSES 

To provide a pictorial representation of a typical GASP 
flight, a case study analysis (adapted from ref, 9) is 
included here. Also, the vertical mean ozone profile and 
latitudinal ozone distribution obtained from analysis of the 
March data are shown. 



Case Study Analysis of Flights on 11-13 March, 1975 

The analyses of the tropopause structure and ozone mixing 
ratios for a series of flights from San Francisco (SFO) to 
Frankfurt, Nest Germany (FRA) via the Far and Biddle East 
are shown in figure 2. The tropopause pressure at each GASP 
data location was obtained by space interpolation from the 
NMC field for the reporting period closest to the flight 
time. Thus, all 0000 GMT tropopause pressure fields are 
assumed valid from 1800 to 0600 GMT, and all 1200 GMT fields 
are assumed valid frco 0600 through 1800 GMT. The structure 
of the tropopause is generally continuous between successive 
analyses for these dates, although a major isobaric 
discontinuity appears between the 0000 GBT and 1200 GBT 
tropopause pressure data on the 13th. This discontinuity is 
evident over the Middle East. 

The GASP ozone and temperature data indicate 
stratospheric flight across the Pacific, the Middle East, 
and South Central Europe. An independent indicator of 
stratospheric flight is provided by the difference between 
the flight altitude and the tropopause height. Here the 
flights across the Pacific and Europe are indicated to be 
stratospheric, in agreement with the GASP data. The 
stratospheric penetration over the Middle East is not 
substantiated by the NBC tropopause height data. This is in 
the region of the discontinuity in the tropopause pressure 
fields, and may be a consequence of a locally varying or 
discontinuous tropopause which the NBC fields would not be 
expected to model correctly. 


For each of the stratospheric flight segments in figure 
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2, ozone naxima are observed as the vlnd direction shifts 
from the southwest quadrant into the northwest quadrant as 
the aircraft flies west. Thus each of the three recorded 
maxima occur within regions of cyclonic wind flow. 


Vertical Ozcne Profile and Latitudinal Distribution 

For three flights around the world (11-13 March, 17-19 
March, and 25-27 March), the variation of the mean ozone 
mixing ratio as a function of the difference between the 
fligUt (pressure) altitude and the NHC-derived tropopause 
pressure is shown in figure 3. Unbiased standard deviations 
are superimposed upon the histogram as well. Here the mean 
values and standard deviations have been obtained for 20 mb 
intervals above and below the tropopause. The increase of 
the ozone nixing ratios beginning at, or near, the 
tropopause is evident as is a coincident naxinun in the 
standard deviation values. 

A zonal display of GASP ozone data above 10 kn for all 
the March flights is shown on figure 4. This indicates the 
high degree of variability of the measurenent s, particularly 
poleward of 30 degrees N. For coaparison. North American 
March mean ozone distribution at 10.0 and 12.5 ka are also 
shown. These curves were calculated fron tabulations of 
molecular concentration given in reference 10. 


CONCLUDING REMARKS 

Atmospheric trace constituents in the upper troposphere 
and lower stratosphere are now being measured as part of the 
NASA Global Atmospheric Sampling Program (GASP), using fully 
automated air sampling systems operating on 747 airliners in 
routine commercial service. Ozone nixing ratio, static air 
temperature, and wind speed and direction data obtained 
during several flights of a GASP-eg uipped Pan American 747 
airliner in March 1975 are now available. The height of the 
tropopause obtained from NMC data archives for the dates of 
the GASP flights are included as a supplement to the GASP 
data. These data may be obtained as GASP tape VL0001 from 
the National Climatic center. Federal Building, Asheville, 
NC, 28801, Plight routes and dates, instrumentation, data 
processing procedures, and tape specifications and formats 
are discussed in this report. 
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TABLE I - GASP PLIGHTS ON TAPE VL0001 


Tape 

Fligh't 

Depa rture 

Data Acgulsition 

Flight 

Nuffibec 

Route 

Date 

Tine (6HT) 

1 

SFO-HHD 

3/11/75 

2013-0618 

2 

HND-HKG 

3/12/75 

0925- 1315 

3 

HKG-BKK 


1512-1702 

4 

BKK-DEL 


1850-2205 

5 

DBL-KHI 


2350-0105 

6 

KHI-BEY 

3/13/75 

0256-0648 

7 

BBY-IST 


0823-0918 

8 

IST-PRA 


1042- 1242 

9 

JPK-SNN 

3/14/75 

0620-1103 

1w 

LHR-JFK 

3/15/75 

1156-1806 

1 1 

JPK-LHR 

3/16/75 

2353-0523 

1? 

LHB-FRA 

3/17/75 

0808-0833 

13 

PRA-IST 


1036-1226 

14 

IST-EEY 


1407-1521 

15 

BEY-KHI 


1659-2006 

16 

KHI-DEL 


220 1-2236 

IV 

DEL-BKK 

3/18/75 

0048-0328 

18 

BKK-HKG 


0537-0717 

19 

HKG-UMD 

3/19/75 

0128-0411 

20 

HHD-SPO 


0640-1440 

21 

SPO-SEA 

3/20/75 

1824-1919 

22 

SEA-LHH 


2143-0529 

23 

LHR-JPK 

3/21/75 

1059-1724 

24 

JPK-LHB 

3/22/75 

1450-2033 

25 

LHR-JPK 

3/23/75 

1040-1345 

26 

JPK-LHR 

3/25/75 

00 36-0607 

27 

LHR-PRA 


0817-0842 

28 

PRA-IST 


10 34- 1229 

29 

IST-BEY 


1410-1455 

30 

BEY-THR 


1711- 1831 

31 

THR-DEL 


2034-2304 

32 

DEL-BKK 

3/26/75 

01 10-0400 

33 

BKK-HKG 


0551-0736 

34 

HKG-HNC 

3/27/75 

0141-0427 

35 

HND-SFO 


0639-1150 

36 

SPO-LAX 


1741- 1801 

37 

LAX-SPO 

3/2 8/75 

1735-1800 

38 

SPO-LAX 

3/29/75 

14 19- 1444 

39 

LAX-GUA 


1709-2039 

40 

GUA-PTY 


2251-0006 

41 

PTY-niQ 

3/30/75 

0240-0354 

42 

MIQ-GUA 


1139-1419 

43 

GUA-LAX 


1656-2031 
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Figure 1. GASP Flight Routes for Tape VLOOOl; Pan Am (N655PA), 11-30 March 1975. Numbers is parentheses 
indicate" number of flights for each route. 
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Variation of ozone mixinq ratio and standard deviation with pressure 
intervals above and below NMC tropopause. 
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LATITUDE, DEG N 

Variation of ozone mixing ratio with latitude. GASP data shown are for 
M^rrh 1975. at fliqht altitudes from 10 to 12 km. 
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APPENDIX A - Specifications for GASP Archive Tapes (VLXXXX) 


GENERAL 

1* Tapes are written in EBCDIC forsat usiig nine track 
tapes. 

2. Tape density is 800 BPI. 

3. Physical records (blocks) are 4096 bytes. 

4. The tapes are unlabeled, with 2 files, a GASP data file 
and a tropopause pressure data file. 

GASP DATA PILE 

1. The GASP data on the tapes is grouped and identified by 

flxghts (takeoff to landing). Each flight begins with a 
logical FLHT record (flight identification data) , which 
is followed by logical DATA records (one for each data 

recording lade during the flight). FLHT and DATA 

records are 512 bytes, hence there are 8 logical records 
per physical record (block) . 

2. A FLHT record will always be the first logical record in 

a block. However, every block need not begin with a 
FLHT record (i.e. , if there are sore than seven DATA 

records in a flight) . If the FLHT record plus the 

available DATA records for a flight do not fill an 

integer nuaber of blocks, the unused logical records in 
the final block are padded with zeros creating PADD 
records. The diagraa below shows how several short 
flights would be blocked. 


Block 12 3 



P D D D D P 

P 

FDD 

D D D D 0 

D D P P P P 

P P 

Logical 

Record 

1 2 3 4 5 6 7 

8 

1 2 3 

4 5 6 7 8 

1 2 3 4 5 6 

7 8 

Block 

4 



5 

6 



FDDDDDDD DDDDDDDD FDDDDDDP 


Logical 

Record 


1 2 3 4 5 6 7 8 


1 2 3 4 5 6 7 8 


1 2 3 4 5 6 7 8 
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where P is a FLHT record 

D is a DATA record 

P is a PADD record 


3* The first four bytes in each logical record identify the 
record type as FLHT, DATA, or PADD. Detailed 

specification of the paraaeters and formats for FLHT and 
DATA records are given in Table A- I and A-II 

respectively. 

a) In each FLHT record, the number of DATA records to 

follow is given by NDATA (Bytes 78-81) , and the 

number of blochs in the flight is given by NBLOCK 

(Bytes 82-84) . 

b) In the last DATA record of each flight, LBFLG 
(Byte 5) = "L"; for the last DATA record on the 
tape, LBFLG = "T*'; for all other DATA records, 
LBFLG =" «. 

Note: DATA records with LBFLG / » ” will be followed by 

PADD records if the physical record (block) is not 
complete. 


TBOPOPAUSE PRESSORE DATA FILE 

1. Following the GASP data, in a separate file, tropopause 
pressure data for the dates of the GASP flights are 
included. Data are currently available from the 
National Meteorological center (NMC) twice daily for 
4225 locations in the Northern Hemisphere. Coordinates 
for these data are the NMC 65X65 square matrix grid, 
transformed from a polar stereographic map of the 
Northern Hemisphere. 

2. Each 65X65 tropopause pressure array is written as seven 
TRPR records. Each TBPB record is a physical record 
(block), and is the sane length as the GASP physical 
records (4096 bytes) . All TRPR records contain 
identification information. Specifications and formats 
for the TRPR records are given in Table A-III. 


Table A-I Foraat for FLHT Records 
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268-273 DFPTA P6.1 Dev/frost point tenperature (deg C) 

274-279 tfVHBA P6. 1 Hater vapor nixing ratio (PPHH) 

280 DPTAGA Al Tag for DPPTA and HVHRA* 

281-286 COAPG F6. 3 Carbon nonoxide data (PPH?) 

287 COTAGA Al Tag for COAVG* 
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the **H” tag is used whenever data is not available 
oc an instrunent is in a calibration node. 
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APPKNDIX B - LATITUDK AND LONGITDDF. FR01 NMC JCC R i. I N AT ES 

Th»^ t ropopfi ur.« jjtossuro tiata lucludod on (;asp taps 1 1 
(jiven odch of the 422S points on the N*1C 6S X hS *jr i .1 , ^ 
square matrix transfotaed from a ^olar stereoyraphic map of 
the Northern Hemisphere. In the hnc coordinates the North 
Pole is the point (33, 3 3), with the 10 deq E - 170 dej V» 
meridian given by the line YJ = 33, and the 100 dey E - HO 
deg W Qieridiin given by the line XI = 33. The 
transf oraatioii from this coordinate system to latitude (dej 
N or S) and longitude (dug E or M) is as follows: 

Let r2 = ((xI-33)2 ♦ (YJ-33)2)/k| (M) 

where Rg = 3 1. 20431 348 


The Latitude ({> (deg) is given by 


0 = (1B0/7I 

) arcsin ( (1 

-r 2)/ (1 iRS) ) 

(A? 

If ()) > 0, 

LAT = <;- 

and LATAG = "N" 


It <}) < 0, 

lat = - <i> 

and LATAG = "S” 



The Longitude 6 (deg) is given by 


0 = “(10 ♦ (IHO/tt ) arctari ( (Y.J- 3 ^)/ (XI- n ) ) ) (A'l) 
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LONG = 
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ind 

LONG'T’AG 
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170 , 
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